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ABSTRACT 

In active galactic nuclei (AGN)-galaxy co-evolution models, AGN winds and outflows are often invoked to explain why super-massive 
black holes and galaxies stop growing efficiently at a certain phase of their lives. They are commonly referred to as the leading actors 
of feedback processes. Evidence of ultra-fast (a > 0.05c) outflows in the innermost regions of AGN has been collected in the past 
decade by sensitive X-ray observations for sizable samples of AGN, mostly at low redshift. Here we present ultra-deep 'KMM-Newton 
and Chandra spectral data of an obscured (Nh~ 2 x lO^^* cm“^), intrinsically luminous (L 2-10 kev~ 4 x lO'^ erg s“*) quasar (named 
PID352) atz~ 1.6 (derived from the X-ray spectral analysis) in the Chandra Deep Field-South. The source is characterized by an iron 
emission and absorption line complex at observed energies of £ » 2 - 3 keV. While the emission line is interpreted as being due to 
neutral iron (consistent with the presence of cold absorption), the absorption feature is due to highly ionized iron transitions (FeXXV, 
FeXXVI) with an outflowing velocity of 0.14/g°gC, as derived from photoionization models. The mass outflow rate — 2 Mq yr“' - is 
similar to the source accretion rate, and the derived mechanical energy rate is ~ 9.5 x 10'*'* erg s“', corresponding to 9% of the source 
bolometric luminosity. PID352 represents one of the few cases where indications of X-ray outflowing gas have been observed at high 
redshift thus far. This wind is powerful enough to provide feedback on the host galaxy. 

Key words. Galaxies: active - Galaxies: nuclei - (Galaxies) quasars: general - X-rays: galaxies 


1. Introduction 


According to some AGN/galaxy co-evolution models, along the 
cosmic history of galaxies there is a dust-enshrouded phase 
associated with rapid SMBH growth and active star forma- 


tion, largely trigge r ed by galaxy mergers a nd encounter s (e.g. . 
Silk & Rees I [T^ loTMatteo et al. I I 200 I iMenci et aTI l200a 


Ho pkins et al. IZubovas & King I I2OI2I : II amastra et al. I 

2OI3I) . At the end of this obscured phase, massive quasar- 
driven outflows blow away most of the cold gas reservoir, cre¬ 
ating a population of “red-and-dead" gas-poor elliptical galax¬ 


ies (e.g., ICattaneo et al. I l2009l) . Support for this picture - at 
least for the most luminous AGN - comes from observations of 


Sturm et al. 11201 1 

; Cicone et al. 112014). as well as of neutral and 

ionized gas (e.g.. 

Nesvadba et al. 1 2008; Alexander et al. 20101: 

Harrison et al. 1 20121 20141 iGenzel et al. 1 120141 Pernaetal. 

2015; Brusa et al! 20151). On smaller scales, ultra-fas 

t out- 

have 

AGN 

flows (UFOs, wi 
been clearly dete 
at low redshift (e 

th velocities typically up to 0.1-0.4c 
cted in X-rays in a sizable sample of 

i.s.. Reeves et al. 1 l2003l iPounds et al. 

2003; 

Tombesi et al.ll201fll l2012al 2014 12015; iGiustini etal. 

2011; 
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Patricket^ 1201 2Hr.offord et al. 1201 .311201 .‘rHKing et al. 

Nardini al. 1120151 : i Ballo et aT] 2015 : see also ICappill2006l for an 

early review, and lFabianll2012l for a more recent one) and also 
in a li mited number of high-redshift quasars (e.g., Qiartaset^U 


llMl2^l2^l2()fl[Saez et al. l2bo^lLanzuisi et al. I2012h . 

Although original claims o f blueshifted absorption features 
were strongly debated (see IVaughan & Uttlevi 2008h . recent 
high photon statistics in X-ray spectra of both radio-quie t 
(e.g., iTombesi et al. I l2()Iol lIoTll mml iGofford et al. 1 12^ 


and r adio-loud AGN (e.g.. iGofford et afl 12013 : Tombesi et al. I 

1201 4l) have undoubtedly shown that these outflows are not rare 
events at X-ray wavelengths, with more than 40% to 50% 
of the investigated Seyfert galaxies in the local Universe be- 
ing associated with such highly energe tic phenomena (e.g., 
iTombesi et al. II2010I: IGoff ord et al. Il2013h . According to mod¬ 
els (e.g., Zubovas & Navakshin 2014ll . fast ionized winds may 
lose most of their kinetic energy after shocking the interstellar 
medium (ISM), thus cooling efficiently and transferring their 
ram pressure (hence momentum flux) to the ISM. In this con¬ 
text, large-scale out flows would be produced by the accelerated 
swep t-up gas (e.g., iKingll^lObt iFaucher-Giguere & Ouatae^ 
1201 ^ . 

In a study based on the X-ray spectral analys is of a sample 
of Type 1 Seyfert galaxies, ITombesi et al. I (1201 3l) have recently 
suggested that X-ray warm absorbers (WA) and UFOs are two 
aspects of the same fundamental physical process, which is re¬ 
lated to wind acceleration likely occurring in the accretion disk 
(see also iFukumura et ^l2014ll . In particular, there is a sort of 
“continuity" in the physical properties of the high- (UFO) and 
low-velocity (WA) winds, with the former being produced in the 
inner regions of the accretion disk, characterized by higher ion¬ 
ization and column density, and the latter associated either with 
the outer parts of the accretion disk or with the AGN torus. While 
WAs have typica l velocities below ~ 1 000 km s~* and low ki¬ 
netic po wer (kg., Bhj^tin_etalj200^ McKernan et al. l2007l but 


see also ICrenshaw & Kraemer| 2012 ). UFOs are characterized 
by velocities > 0.05c up to ~ 0.6c (this h i gh va lue is observed 
in, e g., APM 08279-H52 5.5. LSaez et al. I pOO^ IChartas et al. I 
l2009l and HS 1700-t6416. lLanzuisi et al. f 2012h and consider¬ 
ably higher mechanical energy, so the y possibly exert a signif¬ 
icant impact on the host galaxy (e.g., Hopkins & Elvis J 201ft 


lKing|2010bl:lGaspari et ak 1201 lUZubovas & Navakshin 120141) . 

Regardless of the exact launching site for the wind, the ob¬ 
served outflow properties require an efficient process of ac¬ 
celeration, such as radiation pressure through Compton scat¬ 
tering and, more importantly, mag n eto-hydrodynamical pro¬ 
cesses (e.g., K^_g_^Ppunds_ 200^ ftog^^ICffilmM 20041: 


ISim etal.' 20081: Ohsugaet al.l2009tl^kumura et al. 11201 ol: see 


also Higginbottom et al. Il2014ll. 

From the innermost regions of AGN (probed by X-rays) 
to the outer galactic scales (probed by observations in molec¬ 
ular lines), all of these outflows may provide sign i ficant feed¬ 
back into the q uasar host gala xy (e.g.. lKin^l2010al: iKing et al. I 
1201 It see also iPounds 1 1201 4l for a recent rev iew) and may be 
responsible for quenching star formation (e.g.. lCano-Dfaz et al. I 
l2012tl&esci et aini201.5 ) and, if the momentum is maintained 
through the shocks, for setting up the local black hole mass vs . 
bulge relation (e.g .. lKing Il20101x IZubovas & Navakshin II2014I: 
ICosta et ain l2014ll . However, despite the increasing observa¬ 
tional evidence of outflows at different scales, the link between 
UFOs and large-scale outflows is today basically unknown, es¬ 
pecially at high redshift. In this regard. X-ray emission offers a 
powerful probe of the inner most regions of the AGN where fast 
outflows are launched (e.g.. lProga & Kallman 1120041) . 


In this paper we present the intriguing properties of the 
source XMMCDFSJ033242.4-273815, whose XMM-Newton 
and Chandra spectral data are consistent with the presence of an 
emission-plus-absorption line, although at a different statistical 
significance (lower in Chandra). The emission line is interpreted 
as neutral iron Ka emission, and the absorption line is most 
likely due to a blueshifted highly ionized iron transition asso¬ 
ciated with outflowing gas with a velocity of ~0.1c. The spectral 
“peculiarity" of this source came out in the search for obscured 
AGN at high redshift selected in the XMM-Newton observa¬ 
tions of the Chandra Deep Field-South (XMM-C DFS) by means 
of restframe X-ray colors (llwasawa et al. ll20T^ . This source is 
listed as n umber 352 in the 3M s XMM-Newton 2-10 keV source 
catalog bv iRanalli et al. I (l2013l) : hereafter, we refer to this source 
as PID352. In the following, we present the analysis of these 
iron features based on the whole dataset of XMM-Newton and 
Chandra observations, carried out over a time interval of ~10 
years. Most of the spectral results presented in this paper rely on 
XMM-Aewfon data, while Chandra provides further support to 
these results, as well as an accurate source position. 

Currently, there a re only two photome tric redshift estimates 
for PID352, one from iTavlor etXl (l2(^. ^ (95% 

confidence level), and the other from Hsu et al.l (1201# . Zph = 


1.31 


+ 0.78 

- 0.74 


(99% confidence level). 


2. X-ray data 

The CDF-S was observed with XMM-Newton with 33 exposures 
over the years 2001-2002 and 2008-2010, for a nominal expo¬ 
sure time of ~3.45Ms. Full details on the pointing strategy, data 
cleani ng, and derived catalogs are presented by iRanalli et al. I 
(l2013h . In the analysis presented in this work, we use data taken 
from the three EPIC cameras (pn, MOSl, and MOS2) after ap¬ 
plying a sigma-clipping procedure to the event files in order to 
Alter high-background flaring intervals. Given the large off-axis 
angle of PID352 in the summed XMM-Newton mosaicQ the fi¬ 
nal effective (i.e., after Altering and accounting for vignetting 
effects) exposure time at the source position is ~800 ks (aver¬ 
aged over the three cameras). The resulting 2-10 keV position 
of PI D352 is (ra.dec')^(03:32:42.46.-27:38:15.7: lRmalli et al. I 
l2013h . Spectral data were extracted from the three mosaics of pn, 
MOSl and MOS2 using circular regions of 10" radius (corre¬ 
sponding to C 60% of the encircled energy fraction at the source 
position) in order to maximize the signal-to-noise ratio in the ex¬ 
traction regions. Background was chosen from nearby, source- 
free circular regions of radius 19", 25", and 25" for the pn, 
MOSl, and MOS2, respectively. Once normalized to the source 
extraction regions, the background accounts for ~30-35% of the 
overall counts. The final number of net (background-subtracted) 
source counts in the ~0.3-7 keV band is 2560, 690, and 1250 in 
the pn, MOSl, and MOS2 cameras, respectively. 

PID352 was also observed in both of the Chandra CDF-S 
and Extended-CDF-S (E-CDE-S) exposures. CDE-S data were 
taken in the periods 1999-2000, 2007, and 2010 with 54 point¬ 
ings using the ACIS-I array for a total on-axis exp osure time of 
~4Ms; see lLuo et al. I (1200# and IXue et al. I (1201 Ih for more de¬ 
tails on the observational strategy and d ata reduction. PID 352 
is associated with source XID 571 in the IXue et al. I (1201 ih cat- 

* The exact off-axis angle in the final XMM-Newton mosaic is diffi¬ 
cult to quantify, since the nominal right ascension and declination of 
the individual pointings varied substantially, with the highest difference 
being between the 2001-2002 observations and the most recent 2008- 
2010 observations. 
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alog [with coordinates (ra,dec)=(03:32:42.73,-27:38:15.6)]. As 
a consequence of the large off-axis angle of the source (>10') 
and the different roll angles of the CDF-S exposures, the source 
falls in the CDF-S in only 13 out of the total 54 Chandra point¬ 
ings. The count distribution seems to be slightly displaced (by 
~1.1") with respect to the X-ray source position in the cata¬ 
log, probably because of the off-axis position (hence elongated 
shape) of P1D352 in the 4Ms exposure. Once this offset is taken 
into account, the source is perfectly coincident with a red galaxy 
(see ®. The final vignetting-corrected source exposure time is 
~430ks. The Chandra spectrum was extrac ted from individua l 
pointings using the Acts Extract software (iBroos et al. Il20icil) 
and accounting for the different point spread function (PSF) 
size and shape at the source position in each observation. The 
total number of net source counts in the ~0.5-7 keV band is 
~450, whi c h is fu lly consistent with the number reported in the 
IXue et al~l (1201 Ih catalog. The overall contribution of the back¬ 
ground (chosen from nearby source-free regions) is limited to 
less than 10% of the source plus background counts within the 
source extraction region. We note that at ~8.7" from its position 
there is another fainter (by a factor of ~ 3. 5 in the 0.5-8 keV 
band) X-ray source, listed as XID 568 in IXue et al~l (1201 ih . 
The good agreement in terms of X-ray spectral results from 
'KMM.-Newton and Chandra data suggests that the contamina¬ 
tion by XID 568 to PID352 in XMM-Newton data is most likely 
marginal. 

PID352 was also detected in the E-CDF-S mosaic (four 
pointings with a nominal exposure of ~250ks each; observa¬ 
tions were taken in 2004) and was reported as source 437 
in the iLehmer et al. I (l2005l) catalog. At the source position 
(ra,dec)=(03:32:42.63,-27:38:16.1), which corresponds to an 
off-axis angle of 5.9', the effective exposure time is ~190ks. 
The source spectrum was extracted from a circular region with 
a radius of ~4.5", while the background spectrum was extracted 
from a circular source-free region with radius 19.5". The im¬ 
pact of the background is 4% in the source extraction region. 
The total number of source net counts in the ~0.5-7 keV band 
is ~260, so it is fully consistent with the number reported in the 
iLehmer et al. I (l2005b catalog. 

In the following analyses, XMM-Newton and Chandra (both 
CDE-S and E-CDE-S) results will be presented separately. Over¬ 
all, XMM-Newton data provide good counting statistics, while 
Chandra data offer good source positioning, low background 
contamination, and independent (though statistically limited) 
support for XMM-Newton spectral results. 

X-ray spect ral analysis is carried out using xspec vs. 12.6.1 
(lArnaud II 19961) . Errors are quoted at the 90% confidence level 
for one interesting param eter (i.e., either A^^=2.71 or AC=2.71; 
lAvni 11197^ ICash 1119791 depending on the adopted statistic). All 
spectral fits include absorption due to the line-of-sight Galacti c 
col umn density of Nh == 7 x 1 0^^ cm“^ (iKalberla et al. Il200^ . 
and lAnders & Grevessel (Il989l) abundances are assumed. Here¬ 
after we adopt a cosmology with //o=70km s“' Mpc^*, Qa-0.7, 
and Om=0.3. 


2.1. XMM-Newton spectral results 

The XMM-Newton spectra were binned at a signal-to-noise ra¬ 
tio of at least 3. We have verified that such binning satisfies 
the criterion of the minimum number of photons/bin typically 
adopted to apply statistics (e.g., 20-30 counts per bin). We 
have also checked that the results reported hereafter are recov¬ 
ered (within the errors) once a “standard" binning is adopted. 
In the following, we adopt “phenomenological" models to re- 
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Fig. 1. Plot of the observed 1.3-7 keV band spectrum of PID352 divided 
by the best-fitting powerlaw model (with a flat photon index, see text for 
details). The spectral data were made by combining the XMM-Newton 
EPIC pn, MOSl, MOS2, and Chandra ACIS-I data. The powerlaw fit 
was performed to the data from all these instruments jointly. The most 
relevant spectral features are indicated in the figure. 
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produce the XMM-Newton data, then in 0 we present results 
obtained with a more self-consistent and physically-motivated 
model. 

Eitting the XMM-Newton data with a single powerlaw and 
Galactic absorption produces statistically significant data-to- 
model residuals (;^fVd-0.f=484.5/276, where d.o.f. is the num¬ 
ber of degrees of freedom). Erom the spectral deviations shown 
in Eig. [U (where the assumed model is a powerlaw), it is clear 
that additional spectral complexities are present at observed en¬ 
ergies of ~2-3 keV, namely an emission and an absorption fea¬ 
ture, and an absorption edge. The notably flat photon index, 
F = 0.45 + 0.05, is suggestive of obscuration, and an additional 
component (parameterized by a second powerlaw) is required in 
the soft band. 

Thus we adopted a double powerlaw model with free pho¬ 
ton indices; one of these components is absorbed by cold matter. 
The resulting fit provides ;if^/d.o.f=345.7/273. The two photon 
indices are consistent (within ~ Tier) considering their errors 
(Esoft = 1-94 /;q 34 and Fhard = l-53/o2o)> while the column den¬ 
sity (pha model in xspec) at z-0 is Nh = (1.53/;g3j) x 10^^ cm“^. 
If the soft X-ray emission is interpreted as scattering, this mod¬ 
eling results in a ~6% of scattering fraction, which is a bit high 
but possible based on the observed range in X-ray spectra of ob¬ 
scured AGN (e.g.. iLanzuisi et aini2015l) . Since we are mostly 
interested in the source emission and absorption features and 
the hard X-ray emission, we do not investigate this issue further. 
Then we modified the previous model by leaving the redshift of 
the absorber free to vary. We obtain z = 1.59/^23 and powerlaw 
slopes consistent with the values reported above within their er¬ 
rors; the derived column density is Nh = (2.87/^^^) x 10^^ cm“^. 
The resulting fit quality is^^/d.o.f=301.7/272. The emission and 
absorption features seem to be too separated in energy to be a P- 
Cygni profile. Therefore, to reproduce these features, we start by 
adding a narrow (cr = 10 eV) emission line to the previous mod¬ 
eling, leaving the redshift as a free parameter. Constraints on red¬ 
shift come from the emission line (assumed to be associated with 
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Observed Energy (keV) 



Fig. 2. (Left panel) 'XMM-Newton spectrum of PID352 fitted with a double powerlaw (one absorbed) plus one emission line and one absorption 
line. The black/blue/red lines and datapoints refer to pn/MOSl/MOS2 data. Emission and absorption lines are both marked. In the bottom panel, 
the data-to-model residuals are shown in units of a. (Right panel) Best-fitting 'KMM-Newton model (thick solid line), with indications of the 
spectral components (dotted lines). In particular, we note that the iron edge is at a ditferent energy with respect to that of the absorption feature. 




Fig. 3. XMM-Newton spectral results: (left panel) redshift vs. FeKo- line normalization. For comparison, the redshift estimate obtained by Chandra 
data is shown in the inset; (right panel) observed-frame absorption-line energy vs. normalization. In all panels, contours represent the 68, 90, and 
99% confidence levels for two parameters of interest. 

Table 1. Best-fitting X-ray spectral parameters of PID352 derived from XMM-Newton and Chandra (CDF-S and E-CDF-S). 


Data 

Fsoft 

rhaid 

Nh 

(xl0®cm-2) 


(eV) 

Eabs.line 

(keV) 

(eV) 

F2-IO keV 

(erg cm“^ s“') 

L 2-10 keV 

(erg s“‘) 

xHC)ld.o.f 

X 

c 

1.73+0.31 

1.58!|!| 

1 03+°'® 
l.OJ 0.26 

2 30^°'^* 

^•^^-0.34 

1.59+0.03 

1.62+0.05 

76+35 

120!“ 

2.94+0.04 

9 oc+0.16 

^•OJ_0.07 

-82+“ 

-65^ 

2.6-3.1-2.8 
2.4-3.6 

3.5-1.2-3.8 

1.9-2.8 

275.3/268 

350.0/480 


Notes. X: XMM-Newton (pn, MOSl, and MOS2); C: Chandra (CDF-S 4Ms-hE-CDF-S data). Both line energies and equivalent widths are reported 
in the observed frame, while the column density is given at the reported best-fitting redshift solution. Fluxes are reported in the observed 2-10 keV 
band (units of 10“'^^), while luminosities are intrinsic (i.e., corrected for obscuration) and in the restframe 2-10 keV band (units of 10“^). The 
three fluxes and luminosities reported for XMM-Newton refer to pn, MOSl, and MOS2 data, respectively, while for Chandra the two fluxes and 
luminosities refer to CDF-S and E-CDF-S data, respectively. Both emission and absorption lines were modeled using a Gaussian feature of width 
o-=10eV. 


the neutral restframe 6.4 keV iron transition, which is expected 
given the presence of cold absorption) and the iron Ka edge (at 
the restframe energy of 7.1 keV, seen at an observed energy of 
~ 2.7 keV in Fig.lTJ associated with the absorber. Including the 
line (observed equivalent width EWobs = 76+35 eV, correspond¬ 
ing to TtVPrest ~200 eV at the derived redshift of z = 1.59 + 0.03) 
improves the quality of the spectral fit (;^f^/d.o.f=289.3/271, i.e.. 


A;t'“/Ad.o.f=12.4/l with respect to the previous double powerlaw 
plus absorption model). The line EW appears only marginally 
consistent (within errors) with the EW expected in the case of 
transmission through a neutral shell of gas according to the re¬ 
cent M onte Carlo calculations provided by iMumhv & Yaaoobi 
( 2009h. but seems to agree more with past-wor k predictions (e.g., 
lAwaki et al. II 199 it iLeahv & Creighton1ll993h . 
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The X-ray spectral slopes do not change significantly with 
respect to the previous modeling once errors are taken into 
account. The column density at the newly derived redshift is 
Nh = (2.30+°;™) X 10^3 cm-^. 

The following step consists of modeling the absorption fea¬ 
ture with a narrow (10 eV) Gaussian line. This provides a 
spectral improvement of A;if“/Ad.o.f=14/2 (;^f“/d.o.f=275.3/269). 
Monte Carlo simulations indicate that the probability of by 
chance detecting an absorption line providing an improvement 
to the spectral fitting larger or comparable to the one we obtain 
is ~1% (see ^2.3l for details). Most of the parameters related to 
the X-ray spectrum remain unchanged. The observed-frame en¬ 
ergy and EW of the absorption line are E=2.94 + 0.04 keV and 
- 82/35 eV, respectively. 

The best-fitting XMM-A^ewfo« pn, MOSl, and MOS2 spec¬ 
tra are shown in Eig.|2(left panel). The redshift solution is shown 
in Eig.[3(left panel; z = 1-59/q q 5 according to 90% confidence 
contours), while the observed-frame energy vs. normalization of 
the absorption line is shown in Eig. [2 (right panel). The confi¬ 
dence contours indicate that both the source redshift and the ab¬ 
sorption feature are well constrained in the XMM-A^ewfo« anal¬ 
ysis, which is considerably better than using Chandra data (al¬ 
though consistent in terms of line parameters, see ^2.2b . It is 
worth mentioning that the energy of the neutral iron absorption 
edge (E=2.74 keV) differs clearly from that of the absorption 
line (see Eig.|2] right panel). 

At face value, the derived blueshift of the absorption line 
(most likely ascribed to ionized iron resonant absorption), as¬ 
suming the redshift of z = 1.59, is ~0.13c and ~0.09c in case 
of He-like (EeXXV, 6.7 keV) and H-like (EeXXVI, 6.97 keV) 
iron transitions, respectively. The association with either He- 
or H-like iron transitions is the most likely, according to pre- 
vious, higher statist ics studies of AGN showing outflows (e.g., 
iRisaliti et al. 11200^ : see ®for a complete discussion on this 
issue using a more physically motivated modeling of XMM- 
Newton data. 

We have also verified whether different modeling could ac¬ 
count for the source complex emission; in particular, motivated 
by the observed flat photon index (which may be partly related 
also to the presence of radio emission, see ^3.2b . we checked 
for a r eflection componen t by using the pexmon model within 
xsPEC (iNandra et al. Il20()^ . This model includes Compton re¬ 
flection with self-consistently generated fluorescence emission 
lines from iron and nickel plus the EeKcr Compton shoulder. 
We set the pexmon model to only produce the reflection com¬ 
ponent by making the relative reflection parameter R negative 
and by tying its photon index and normalization to the cor¬ 
responding parameters of the primary powerlaw. We fixed the 
high-energy cutoff of the illuminating powerlaw at 100 keV, the 
abundances of the distant reflector to the solar values, and the 
inclination angle to 60 degrees. This model does not provide 
a good fit to the data (;^f^/d.o.f=323.8/270), because the reflec- 
tion/powerlaw normalization ratio is unphysically high and ba¬ 
sically unconstrained, as the entire spectrum of PID352 were 
reflection-dominated (which is also not consistent with the ob¬ 
served emissi on-line EW). Similar results h ave been obtained 
using PEXRAV (iMagdziarz & Zdziarski I [T995h plus an emission 
line. 

The observed-frame 2-10 keV source flux is ~ (2.6 - 3.1 - 
2.8) X 10“'"^ ergcm^^ s“' (where the three values refer to pn, 
MOSl, and MOS2, respectively, owing to their slightly differ¬ 
ent normalizations). Once corrected for intrinsic absorption, the 
restframe 2-10 keV source luminosity is ~ (3.5 - 4.2 - 3.8) x 


lO'*"^ erg s '. The best-fitting 'KMM.-Newton results are summa¬ 
rized in Table [1] 

2.2. Chandra spectral results: CDF-S and E-CDF-S 

We then checked for indications of the iron line complexity in 
PID352 using the available Chandra spectra, well aware that 
the limited photon statistics will prevent us from drawing firm 
conclusions on this issue from Chandra data alone. At first, the 
CDE-S and E-CDE-S spectral data of PID352 were fitted sepa¬ 
rately in order to appreciate, from the admittedly low-counting 
statistics spectra, any possible significant difference, mostly in 
the continuum emission, which might preclude a joint analysis. 
A powerlaw model was adopted at this stage. The datasets pro¬ 
vided consistent spectral results in terms of photon index (within 
errors), while the source flux is a factor ~1.4 higher in the E- 
CDE-S data than in CDE-S data. In all of the following analy¬ 
ses, we proceed with simultaneous X-ray spectral fitting of both 
datasets, leaving the normalizations free t o vary and u sing the 
Cash statistic (i.e., using unbinned data; see lCash Ill979h . 

The flat photon index obtained adopting a powerlaw model, 
E = 0.2 + 0 . 1 , is suggestive of the presence of obscuration 
toward the source, similar to 'KMM.-Newton data analysis. To 
properly investigate this issue, we included an obscuring screen 
in the spectral modeling, which returns a column density of 
( 9 . 6/3 5 ) X 10^' cm“^ at redshift z-0. As expected, the photon 
index becomes steeper (E = 1.2 + 0.3), although not as steep 
as ty pically found in unobs cured AGN and quasars (E a; 1.8, 
e.g., iPiconcelli et al. I l2005h . The shape of the data-to-model 
residuals, especially in the higher statistics CDE-S spectrum, 
are suggestive of some spectral features at observed energies of 
~ 2 - 3 keV, although the relevance of such features is clearly 
limited by the paucity of photons. 

Then we included an emission line, whose best-fitting 
observed-frame energy is 2 . 44/553 keV, consistent with iron Ka 
emission if the redshift derived from XMM-Newton data analy¬ 
sis is assumed. The improvement with respect to the previous fit 
in terms of AC is limited (10 for two less d.o.f., i.e., the line en¬ 
ergy and its normalization, since the line width has been frozen 
to (T=10 eV). 

We note that other data-to-model residuals may reflect the 
presence of an absorption line, similar to what has been ob¬ 
served in XMM-Newton data. Once modeled as a narrow (10 eV) 
Gaussian line, the fit improves slightly (AC/Ad.o.f.=s5/2). The 
observed-frame energy of this additional line is 2. 85/° ™ keV. 

As for the XMM-Newton spectral analysis ( ^2.11) . we then 
decided to model the data assuming that the iron Ka line is neu¬ 
tral (6.4 keV restframe) and leaving its redshift free to vary. The 
redshift solution is therefore determined by the iron emission 
line and the iron absorption edge due to neutral matter. The de¬ 
rived redshift is z = 1.62 + 0.05 (z = 1. 62/553 if one considers 
the 90% confidence level contours plotted in the inset of Eig. [2 
left panel). The column density at the newly determined source 
redshift is Nh = (I.OI/q^) x 10 ^^ cm“^, while the photon in¬ 
dex apparently remains flat (E = 1.0 + 0.3), possibly suggesting 
additional absorption complexity. Assuming E=1.8, the column 
density becomes ~ 2 x 10“^ cm“^. The observed-frame emission 
(absorption) line EW is 120 eV (-65) eV, with large uncertain¬ 
ties (see Table[T]i owing to the low statistical significance of both 
lines. The best-fitting X-ray spectrum (rebinned at the 3cr level 
for presentation purposes) is shown in Eig. 0] where E-CDE-S 
data are shown in red and CDE-S data in black, both normalized 
by the effective area. The energy of both features is marked. At 
the redshift of z = 1.62 and assuming that the absorption fea- 
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Observed-frame Energy (keV) 

Fig. 4. (Top panel) Chandra CDF-S and E-CDF-S spectra of PID352, 
normalized by the effective area and fitted with an absorbed power- 
law, plus one emission line and one absorption line. The red datapoints 
refer to the E-CDF-S data, and the black ones to the higher statistics 
CDF-S data. The best-fitting model is shown as continuous curves. The 
observed-frame energy of the emission and absorption lines are also 
shown. All of the data have been rebinned to the 3cr level for presenta¬ 
tion purposes. (Bottom panel) Data - model residuals once both datasets 
are fitted with an absorbed powerlaw without including any feature. 


ture is due to either He-like or H-like iron transitions, we derive 
blueshift velocities of ~0.1 Ic and ~0.07c, respectively. 

The flux of the source in the observed 2-10 keV band is ~ 
(2.4 - 3.6) X 10“'^ erg cm“^ s“', corresponding to a restframe 
2-10 keV intrinsic (i.e., corrected for the measured absorption) 
luminosity of ~ (1.9 - 2.8) x 10^“* erg s“' (where the two values 
refer to CDF-S and E-CDF-S data, respectively). The best-fitting 
Chandra results are summarized in Table [1] 

2.3. What XMM-Newton and Chandra tell US about PID352 

From the X-ray spectral analyses reported above and summa¬ 
rized in Table [1] it is clear that the XMM-Newton and Chandra 
results are consistent within errors, both in terms of continuum 
emission and emission/absorption features. The XMM-Newton 
flux lies between that of CDF-S and E-CDF-S data. The only sig¬ 
nificant spectral difference is the presence of an additional soft 
component, parameterized by a powerlaw, in the XMM-Newton 
data; it seems plausible that we are able to observe such a com¬ 
ponent in XMM-A^ewfon because of its higher effective area with 
respect to Chandra. The hard powerlaw is flatter in Chandra than 
in XMM-Newton and may explain the apparently higher column 
density derived in the XMM-A^ewfo« spectral analysis (see Ta¬ 
ble [TJ. Once F = 1.8 is assumed, the column densities derived 
from Chandra and XMM-Newton data become similar. We note 
that applying the best-fitting XMM-Newton model to the Chan¬ 
dra data (allowing the normalization of the continuum to be dif¬ 
ferent) provides an acceptable fit. Also the features associated 
with iron have similar properties in XMM-Newton and Chandra 
spectra. 

Assuming the best-fitting results from XMM-Newton analy¬ 
sis (owing to the higher quality of the data) and z - 1.59 + 0.03, 
we are able to derive the following outflow velocities for the gas: 
Wout = 0.13 ± 0.02c and 0.09!° “c in case of He-like (FeXXV, 
6.7 keV) and H-like (FeXXVI, 6.97 keV) iron transitions, re¬ 
spectively. We postpone the discussion on the ionization state of 
the line to Sect. 01 


Although the absorption feature is observed in both XMM- 
Newton and Chandra datasets, one may argue about its statis¬ 
tical significance. To this point, extensive Monte Carlo simu¬ 
lations were carried out by producing 10 000 spectra for each 
EPIC instrument ( pn, MOSl, and MOS2 ) using the fakeit rou¬ 
tine in xsPEC (e.g.. iLanzuisi et aT~ll2013h . Eor each camera, the 
simulated spectra have the same background flux and exposure 
as observed. As input model for our simulations, we assumed the 
best-fitting model obtained in Sect. 12. II deprived of the absorp¬ 
tion feature. Therefore, the significance of the absorption line 
is estimated by computing how many times this feature is re¬ 
covered by chance and produces an improvement in the spectral 
fitting larger or comparable to the one actually measured using 
real data. In practice, we first fitted the simulated spectra with 
the same model used for the simulation with free parameters and 
then introduced an absorption line. The line energy was left free 
to vary in the full 0.3-7 keV range and constrained to be the same 
in all three (pn, MOSl, and MOS2) spectra. To avoid the fitted 
energy of the absorption line clustering around the initial value, 
for each simulation we repeated the fit with an additional absorp¬ 
tion line 68 times with initial energies in steps of 0.1 keV be¬ 
tween 0.3 and 7 keV, and selected the best fit among them. Only 
101 of the 10 000 simulated spectra show a spectral improve¬ 
ment larger than or equal to what is observed {x~ ^ 14). This 
translates into a probability of detecting an absorption line when 
it does not exist of 1.01%. We then used a Bayesian procedure 
to find the probability distribution of the fraction of spurious de- 
tections (see example in pages 20—22 of lWall & .lenkins]l2003h 
from 101 spurious detections out of 10 000 trials. This probabil¬ 
ity distribution was used to estimate confidence intervals in /, by 
determining the narrowest interval around the mode / = 0.0101 
(S. Andreon, private communication) that included the corre¬ 
sponding probability. We obtained / < 0.0132 at 99.73% confi¬ 
dence. 

As IVaughan & Uttlevi (l2008h point out, relatively strong 
lines with large uncertainties (low EW/error ra tio) may suffer 
from the so-called publication biafl (but se e alsolT ombesi et al. I 
l2010h . In the XMM-CDES source catalog (iRanalli et al. 112013h . 
there are only ten sources with more than 5000 net counts, 
PID352 being exactly in the tenth position in terms of number 
of counts. X-ray spectral analysis o f this high-statistics sa mple 
(Comastri et al., in preparation; see llwasawa et al.1l2015l for a 
detailed study of the two brightest AGN) shows no indication of 
absorption features as the one detected in PID352. We can safely 
assume that below this number of counts, the detection of ab¬ 
sorption features would be challenging. Considering the size of 
this XMM-Newton highly reliable “spectroscopic" sample an d 
the discussion reported in §4.4 of IVaughan & Uttlev I (l2008l) . 
we can argue that we might expect 0.1 such detections above 
the ~2.6cr significance of the line. A comparable number is ob¬ 
tained if we consider also the four sources with more than 4000 
net counts. The presence of such a feature in the different EPIC 
cameras and, possibly, in the Chandra data makes us confident 
of the reliability of the line, although its statistical significance 
is admittedly limited. 

Eurthermore, we note that at the energy of the emission and 
absorption features, there are no calibration issues in the effec¬ 
tive area of XMM-Newton instruments. The accuracy is better 


^ The publication bias consists in the fact that only the observations 
with detected features have been reported in the literature. In other 
words, the reported detections of emission/absorption lines may be just 
the “tip of the iceberg", i.e., the strongest or most significant ones from 
a distribution of random fluctuations. 
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than 3% and 2% in MOS and pn cameras, respectively (see 
XMM-SOC-CAL-TN-0018). The main response artifacts (i.e., 
the An M and Ir K edges in the effective area) can therefore not 
be responsible for the observed features. 

3. Multiwavelength data 

3.1. Spectral energy distribution of PID352 

Using the available multiwavelen gth data for P1D352 (from 
MUSYC, ICardamone et al. I l2010l Spitzer, Herschel SPIRE, 
SCUBA, LABOCA, dxvAAzTEC), we were able to characterize the 
source more properly. In part icular, P1D352 is as sociated with an 
extremely red object (ERO. iHston et al71ll988h . because it has 
an R-K color (AB mag) of ~3.8 (5.4 in Vega mag). This class 
of sources represents one of the major c omponents of the s tellar 
mass build-up at redshifts z ~ 1-2 (e.g.. lDaddi et al. l2000l) . The 
source flux density in the observed MUSYC 7 band is a factor ~3 
higher than in the z band, where P1D352 is barely visible (see 
Fig.B the source position is shown as X-ray contours provided 
by Chandra CDF-S data). To have a more comprehensive pic¬ 
ture of PID352 and derive an estimate of the source bolometric 
luminosity, we modeled the source UV to millimeter emission in 
terms of galaxy plus AGN components (see Fig.|6]l. 

The spectral energy distribution (SED) of P1D352 was 
modeled u s ing th e SED-fltting code orig inally developed by 
iFritz et alTI (l2006l) and recently updated by iFeltre et aTI (l2012l) . 
and assuming z - 1 .60 as derived from our X-ray spectral analy¬ 
sis. We note that, within the errors, the X-ray-derived redshift 
is consistent with both current pho tometric-redshift estimates 
dTavlor et al. 11200^iHsu et al. Il2014li . and it is currently the one 
with the lowest uncertainties. The code accounts for both the 
stel lar and the AGN components (for applications of this code , 
see IVignali et al. I l2009t (Pozzi et ahl 1201 2t iGilli et ahl l2014l) . 
The stellar component is composed of a set of simple stellar pop¬ 
ulation spectra of solar metallicity and ages up to ~3 Gyr, i.e., the 
time elapsed between Zform=6 (the redshift assumed here for the 
stars to form) and z=1.6 (the source redshift). This component 
is extincted according to a lCalzetti et al. I(l2000l) attenuation law. 
For the AGN component, we used an extende d grid of “smooth " 
torus models with a “flared disc" geometry (iFritz et al. II2006I) . 
which provides a good description of the AGN mid-lR SED in 
the ca se of sparse photometr ic datapoints (see also lVignali et al. I 
IMil an d lFeltre et al. 112011 . 

The galaxy provides the dominant contribution at short 
wavelengths (dotted line in Fig. |6ll, where the best-fitting torus 
solution is suggestive of an obscured AGN (dashed line in 
Fig. ID, i.e., an AGN providing a limited contribution to the op¬ 
tical emission because of extinction. The amount of obscuration 
toward the AGN derived from the SED fitting is fully consistent 
with the absorption found from X-ray spectral analysis (assum¬ 
ing a Galactic dust-to-gas ratio conversion). However, the prop¬ 
erties of the AGN component that may be derived from the SED 
fitting are not well constrained because of the lack of additional 
mid-IR data besides the Spitzer MIPS 24/rm datapoint. 

This said, we used t he restframe 12.3pm flux density and 
the iGandhi et al. I (l2009l) relation to estimate the predicted in¬ 
trinsic 2-10 keV luminosity of the source. We obtained a value 
of ~ 2.5 X 10^"^ erg s“\ which is ~40% lower than actually mea¬ 
sured from XMM-Newton spectral analysis but overall consis¬ 
tent if one considers the scatter in this mid-lR vs. X-ray correla¬ 
tion. 

Although the sou rce is apparently de tected at 70pm in 
the Herschel catalog (iMagnelli et al. Il2013ll . the presence of a 



Rest —frame wavelength (/xm) 

Fig. 6. Restframe spectral energy distribution for PID352. Available 
datapoints and upper limits (Scr) are plotted as red filled circles and 
downward-pointing arrows, respectively. The solid black line is the 
summed contribution of the AGN (blue dashed) and galaxy (red dot¬ 
ted) components. A redshift of z = 1.6 has been assumed. 


nearby far-IR bright source and the lack of any further detec¬ 
tion of P1D352 at longer wavelengths (including LABOCA and 
AzTEC data at 870pm and 1100pm) suggest to consider P1D352 
as undetected in the far-lR/mm (as shown in the SED plotted in 
Fig.lD. 


3.2. Radio data 

P1D352 has a long history of observations in the radio band. 
From the analysis of VIA 20c m (1.4GHz) and 6cm (4 .8GHz) 
maps, iMiller et al. I (l2008ll and iKellermann et al. I (l2008h report 
the presence of an exten ded radio galaxy, classi fied as a dou¬ 
ble Fanaroff-Riley (FRll, iFanaroff & Rilev II 19741) . whose “cen¬ 
troid" corresponds to the red galaxy associated with P1D352 (see 
Fig.EJ. However, no radio emission is detected at the position 
of the core of this FRII galaxy, although in the 6cm map there 
are hints of emission at that position (~ 0.64 + 0 .37 mJy). In 
the lat est published version of the 20cm catalog bv iMiller et al. I 
(l2013l) . characterized by an average rms of ~ 7.4pJy and a beam 
size FWHMss 2.8" x 1.6", the western and eastern radio lobes 
have integrated flux densities of 46.1 mJy and 36.2 mJy, respec¬ 
tively (with errors of ~ 30 pJy). The radio power at 1.4GHz is 
thus ~ 9.3 X 10^® W/Hz (iBonzini et al. 112012ll2013ll. If we adopt 
the me asurements of the two radio lobes from iKellermann et al. I 
(l2008l) at 6cm, we derive an energy indejfl of a 0.9 - 1.0, 
which is consistent with the relatively old electron populations 
typically expected in radio lobes. 

The source is also reported i n the ATL AS survey (carried 
out with ATCA) bv iNorris et al. I (l2006l) and iHales et al. I (l2014l) 
at 20cm with a different resolution (~ 11" x 15" and ~ 12" x 6", 
respectively). The flux densities reported in the two ATCA-based 
papers for the western lobe are 44.9+9.0 mJy and 55.4+2.8 mJy, 
respectively, while those for eastern lobe are 27.4 + 5.5 mJy and 
38.8 +1.9 mJy, respectively. Not surprisingl y, PID352 is unde- 
tected in the VLBA observations at 20cm by iMiddelberg et al. I 


^ Here the flux density is reported as Fy oc v 
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Fig. 5. MUSYC z-band (left panel) and /-band (right panel) smoothed images with the Chandra contours of source PID352 overlaid. The X-ray 
source is associated with a red galaxy clearly appearing in filters redward of the z band. Each image is ~ 22" by side. North is up, and east to the 
left. 



Fig. 7. CDF-S 0.5-7 keV unsmoothed im age (40" by side) with the 
20 cm contours (from IMiller et al.1l2013h overlaid. X-ray counts are 
most likely ascribed to the core of t he radio galaxy, w here no appre¬ 
ciable radio emission is detected (see lMiller et al. Il2013h . The Chandra 
position of PID352 (see ^for details) is shown as a circle with a radius 
of 2". North is up, and east to the left. 


4. Properties of the outflow: location, mass-outflow 
rate vs. accretion rate, mechanicai energy, and 
efficiency 

After having characterized the X-ray data of PID352 us¬ 
ing a “phenomenological" model, we modeled the XMM- 
Newton data using grids from the photoionizatio n code xstar 
dBautista & Kallniam 120011 : iKallman et al. 1120041) to derive the 
physical parameters for the X-ray outflowing gas (primarily 
its location and the mass outflow rate), which are clearly af¬ 
fected by large uncertainties given the limited significance of 
the absorption line. The relatively large restframe EW of the 
iron absorption feature in the XMM-Newtoii best-fitting model 
(~ 210 eV assuming z = 1.59, see Table[T]), combined with the 
cu rve of growth for hig hly ionized iron transitions (see Fig. 1 
of iTombesi et al. 1 1201 ll) . suggests the need for high turbulent 
velocities. For this reason, we chose, among the xstar grids. 


the one with a turbulent velocity of 5000 km s“*. This value is 
not at odds with the line width (cr), whose 90% upper limit is 
~ 150 eV (equivalent to ~6000 km s“') in XMM-Newton spec¬ 
tra. The free parameters of the xstar model are the column den¬ 
sity (A^h) of the ionized gas, the ionization parameter (defined as 
^ = Tion/n where Lion is the ionizing radiation in the 13.6 eV- 
13.6 keV energy range and nr the measured column density Ah), 
and the observed redshift of the absorber Zo, which is related to 
the intrinsic absorber redshift (i.e., in the source restframe) Zs. 
as (1 H- Zo) = (1 -H Zfl) (1 + Zc), where Zc=l-59. The velocity 
can then be determined using the relativistic Doppler formula 
^ + Za - ((1 -i8)/(l + j0))°'^, where P - u/c is requi red here to 
be positive for an outflow (see iTombesi et al. Il2()lll for further 
details). 


Using xstar, we were able to obtain a relatively good fit for 
our data (;t'^/d.o.f.=273.5/267), and the derived spectral parame¬ 
ters related to the powerlaws, cold absorption, and iron emission 
line are broadly consistent, within their uncertainties, with those 
presented in ^2.11 The resulting column density of the ionized 
gas and its ionization parameter are NH=(2.23/;j gg) x 10“^ cm“^ 
and log(^) = 2.94 /q 28 erg cm s“\ respectively. These values are 
in go od agreement with tho s e found in AGN with X-ray outflows 
(e.g. ITombesi et ani2010t iGofford et al. 1 1201^ ICh^tas et al. I 
120141) . From the redshift of the absorber, zo = 1.24/° we are 
able to estimate Vout = 0.14/° °gC, i.e., the absorption line is ei¬ 
ther associated with the FeXXV Hea transition at 6.70 keV or 
the FeXXVI Fyo' transition at 6.97 keV (or a combination of the 
two). Given the errors associated with the velocity of the out¬ 
flow and the other spectral parameters, it is difficult to clearly 
determine what the exact ionization state of the line is. We note, 
however, that the 6.7 keV solution generally allows for higher 
EW valu es (i.e., more consiste nt with ours) at a given log(^) (see 
Fig. 2 of iTombesi et al. ir201 ih . 


In the following cal culati ons, we adopt the “pres c ription s" of 
ITombesi et al. and ICrenshaw & Kraemer (12012^ see 

also Krongold et al. I (l2007l) and Pounds & Reeves I (l2009l) for 
additional discussion of outflows and related geometry. 


The minimum distance of the outflowing gas can be esti¬ 
mated from the radius at which the observed velocity corre¬ 
sponds to the escape velocity, i.e., - '2-GMshIvI^^i, where 
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\out is the measured outflow velocityQ. To this aim, we esti¬ 
mate the mass of the black hole directly from the SED fitting 
( ^3.11), where the ste llar mass (M* ^ 4.9 x lO" Mq, assum¬ 
ing a ISalDeterlll959l initial mass function) has been converted 
into a black hole m ass adopting the relation reported in Eq. 8 of 
ISani et akl (1201 Ih . The rms in this relation implies an error on 
the derived mass of a factor of ~ 2.4. Erom the black hole mass 
of ~ 5.0 X 10^ Mq, we derive ~ 7.6 x 10*^ cm, which is 
equivalent to ~100 gravitational radii (rg = GMbh/c^)- Consid¬ 
ering the uncertainties on the black hole mass, the range for rmin 
is ( 2.6 - 26) x 10 ^^ cm. 

Using the definition of the ionization parameter, we can 
also estimate the maximum distance of the gas from the cen¬ 
tral source as rmax < Tion/^ Nh of ~ 2 x 10 "^ gravitational 
radii0. At this point, we can provide a crude estimate of the 
mass outflow rate using the formula reported in Sect. 3.1 of 
ICrenshaw & Kraemerl (I2012h : 

Mout = 4m«p/iAHyout^Cg, (1) 

where m-p is the mass of the proton, // the mean atomic mass per 
proton (=1.4 for solar abundances), Ah the column density of 
the ionized gas, Vout is assumed to correspond to the line of sight 
outflow velocity, r is the absorber’s radial location, and Cg the 
global covering factor . Using a sample of local Seyfert galaxies, 
iTombesi et al. I (1201 Ol) estimate the co veri ng factor statistically to 
be ~0 .5 (see also Gofford et al. l201^ and lCrenshaw & Kraemerl 
l 2012 h . In the biconical approximation of the wind, this corre¬ 
sponds to Cg = Q/ 47 r = 0.5, where Q is the solid angle sub¬ 
tended, on average, by the wind with respect to the X-ray source. 
In the following, we assume r as the minimum distance of the 
outflowing gas (r^in), so the derived parameters are to be con¬ 
sidered as lower limits. We note, however, that our object is at 
least an order of magnitude more luminous than the samples of 
local Seyfert galaxies from which Cg was estimated statistically. 
Although we do not know this value exactly at high redshift and 
high luminosities, we can assume that it should be of this order, 
since it is most likely associated with a wide-angle accretion disk 
wind. 

We estimate a mass outflow rate of Mout ~ 1.7 M© yr"'. The 
large uncertainties in the quantities used to estimate Mout imply 
a range for the mass outflow rate of ~ (0.3 - 6 . 8 ) Mq yr“' (Icr 
significance level). The uncertainty in the outflow rate is already 
quite large, and thus all quantities reported below will have a 
similarly large uncertainty. 

As a consistency check, we computed the source accretion 
rate. Using the bolometric luminosity of PID352 (~ 10^® erg s“') 
derived from the SED fitting and ascribed to accretion processes 
(i.e., the mid- IR emission as due t o thermally reprocessed AGN 
emission; see IVignali et al. Il 201 ll for details) and assuming an 
energy conversion factor 77 = 0.1, we obtain Macc ~ 1.7 Mq yr“', 
i.e., of the same order of magnitude as the mass outflow rate. 

At this stage, one might consider whether the outflow me¬ 
chanical energy is sufficient to influence the st ellar p r ocesses 
in t he host galaxy, following the reasoning of iKing I (1201 Oah 
and [Pounds I (l2014t) . Assuming an outflow velocity of 0.14c, 
we obtain a mechanical energy rate E„,ech - ( 1 / 2 ) Mout v^ut ~ 
9.5 X lO"^'^ erg s“' (uncertain by at least one order of magnitude), 
which is only 1.5% of the Eddington luminosity of PID352 
(~ 6.5 X 10"^® erg s“*). Using a Mbuige ~ 500 x Mbh value and 


a dispersion velocity cr » 280 km s“' (as expect ed from the M- 
<T rela tion given the mass of PID352 black hole: lGultekin et akl 
l2009h . the binding energy of the bulge gas, Ebind ~ IWbuigeO"^, 
is ~ 3.8 X 10^® erg. In a Salpeter time, the mechanical energy 
injected by the wind into the surrounding galaxy medium is 
~ 1.3 X 10®°77niech efg, where T/mech is the unknown fraction of 
the wind energy actually transferred to the bulge gas. Although, 
at face value, the energy injected by the outflow could be even 
greater than the binding energy of the bulge, we note that main¬ 
taining such a wind on bulge scale for a Salpeter time appears 
far from being an easy process. We also note that the efficiency 
of the outflow in PID352, (1/2) Mout f^lutl^hoi ~ 0.09, is consis¬ 
tent with the v alues required for efficient AGN feedback drive n 
by winds (e.g., iDi Matteo et al. Il2005l: iHopkins & Elvis II 2 OIOI 1 . 
As such, the wind in PID352 is likely to have a significant ef¬ 
fect on the host galaxy. The relevance of the current result re¬ 
lies on the high redshift of PID352, close to the observed peak 
of the black hole accretion rate density (e.g., iDelvecchio et al. I 
l2014h . Compared to other high-redshift (z - 1.5 - 3.9) quasars re¬ 
ferred to in the literature, the efficiency of the outflow in PID352 
is not so extreme (efficiencies of ~ 0.1, 0.3, and 1.0 were de¬ 
rived for the lensed quasars HS 0810-H2554, PG 1115-1-080 and 
APM 08279H-5255, respectivelv: FChartas et al. II200^|20T1 . but 
our estimate is based on lower-quality X-ray data and suffers 
from large uncertainties. We note that high values for this effi¬ 
ciency may favor mag netic driving as the accele ration mecha¬ 
nism for the wind (e.g.- lBlandford & Pavne Ill982h . We also note 
that the derived wind momentum. Pout - Mout I’out, is consistent 
with Lhodc. 


Prom recent literature (e.g., iTombesi et al. I 120101: 
iGofford et al. 1 1201 3h . the median velocity derived from the 
highly ionized iron transitions observed in local AGN is ~0.1c; 
a similar velocity has also b een observed in h ighly ionized 
oxygen transitions in Ark 564 (iGupta et al. 1120131) . This average 
value is close to the velocity derived for the outflow observed in 
PID352. Compared with other AGN showing UPOs, the pecu¬ 
liarity of PID352 is that it represents one of the few cases where 
an X-ray outflowing wind has been clearly detected at z > 1 
(see a lso IChartas et al. 11200^ l2007l l2009l 12014t iLanzuisi et al. I 
l2012h . Purthermore, to our knowledge, PID352 is not a lensed 
quasar, while most of the z > 1 outflow detections are related to 
lensed systems. Compared to the lensed high-redshift quasars 
with outflows HS 0810-H 2554 (z=1.51: l^rtas et al. I Iml 
and H1413-I-117 (z=2.56: IChartas et al. Il200'7 ). the bolometric 
luminosity of PID352 is a factor of ~3 and ~ (5 - 10) higher^ 
The relatively good constraints for the outflow parameters of 
PID352 are mostly due to the extensive observational X-ray 
campaign carried out in the CDP-S over more than a decade (in 
this regard, see also the claimed but never confirme d outflows 
reported bv lWang et al. ll200^IZheng & Wang Il2008h and by the 
large collecting area of XMM-Aewfo«. 

As a final remark, we briefly discuss the presence of the 
wind in PID352 and its possible link to radio jets, since the 
so urce is classified as an PRII galaxy (see ^3.2b . As discussed 
by ITombesi et al. I (l2014l) . where a minimum wind covering frac¬ 
tion of C « (0.3 - 0.7) was derived, UPOs in radio-loud AGN 
are likely to be covering a significant portion of the sky as seen 
by the central source. Since derived statistically for a sample of 
26 radio-loud AGN (using also the information from the radio jet 


We assume that the gas is ejected at the observed velocity, without 
accounting for possible additional acceleration of the gas in the flow. 

^ The implicit assumption in the r^ax estimate is that the ionizing 
source is seen as pointlike by the absorber. 


^ We used a bolometric correction of 25 to convert the 2-10 keV lumi¬ 
nosity into a bolometric luminosity, and magnification factors of 94 and 
(20-40) to have intrinsic values for HS 0810+2554 and H1413+117, 
respectively. 
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inclination angle to visualize the geometry of the system, assum¬ 
ing that the jet is perpendicular to the accretion disk), this result 
implies that the absorbing material is not highly collimated as 
expected in the case of jets; as a result, an accretion-disk wind is 
favored. Furthermore, in the specific case of another FRII radio 
galaxy, 3C 111, the UFO mass outflow rate (~ 0.1 - 1 Mq yr“', 
i.e., similar to that of PID352) is much greater than the out¬ 
flow funnelled into the jet {Moutj = 0.0005 - 0.005 Mq yr“*; 
iTombesi et al. I I^12bll . implying that the UFO is much more 
massive than the jet, although its kinetic power can be one order 
of magnitude lower. Placed in a broader context, these results 
are suggestive of considerable feedback on to the host galaxy 
by winds and jets (see also the discussion in iMiller et al. 11200^ 
about the co-existence of outflows and jets in radio-loud quasars 
and their connection with X-ray binaries). 

5. Summary 

In this work we have investigated the properties of PID352, a lu¬ 
minous (L 2-10 kev~ 4x10"^"^ erg/s), obscured (Nh~ 2x10 “^ cm“^) 
quasar at zx ~ 1-6 detected in the CDF-S with both XMM- 
Newton and Chandra and associated with the radio-undetected 
core of an FRII galaxy at 20cm, which is optically classified 
as an extremely red object. X-ray spectral data show the pres¬ 
ence of an iron-line complex consisting of an emission plus an 
absorption line. While the former feature (used, in conjunction 
with the iron absorption edge, to estimate the source redshift) 
is probably associated with neutral iron Ku emission, the latter, 
according to a fit with xstar, is ascribed either to the FeXXV 
or to the FeXXVI transition (or a combination of the two) in an 
outflowing gas, implying iiout = 0-14j;QQgC. The outflow veloc- 
ity is within the observed range for local Seyfert galaxies (e.g., 
iTombesi et al. ll^oT^lGofford et al. Il2013h . Monte Carlo simula¬ 
tions indicate that the probability of detecting an absorption line 
when it does not exist and providing an improvement to the spec¬ 
tral fitting larger or comparable to the one we obtain is ~1%. A 
basic and rather uncertain calculation places the minimum dis¬ 
tance of the outflowing gas at ~ 100 gravitational radii from the 
black hole, i.e., on scales typical of the accretion disk. In the 
panorama of AGN with high-velocity outflowing gas, the pecu¬ 
liarity of PID352 consists in being one of the few sources with a 
detected UFO at z > 1 without being a lensed quasar. 

For PID352 a mass outflow rate of ~2 Mq yr“^ (range 0.3- 
6.8 Mq yr^*) is derived, which is consistent with the source ac¬ 
cretion rate. We have also estimated that the mechanical power 
injected by the wind onto the host galaxy may be very high and 
potentially able to disrupt the bulge, provided that a large amount 
of the wind energy is efficiently transferred to the bulge gas on 
long (> Salpeter) timescales, which looks like an extreme sce¬ 
nario. In any case, the wind is likely to have a significant impact 
on the host galaxy. 

An interesting step forward in the study of PID352 can 
consist of detecting large-scale outflows at longer wavelengths. 
This would allow us to investigate the physical link between 
the small-scale high- velocity outflow prob ed by X-rays and the 
outer-scale outflows dTombesi et al. 11201^ . In particular, given 
the source magnitude in the near-IR (AB magnitude ~21.8 and 
21.3 in the J and H bands, respectively), a moderate-length 
(~1 hr) observation at ESO with VLT/XShooter (allowing for 
sensitive, medium-resolution spectroscopy up to the Xband) will 
be sufficient to detect the main restframe optical emission lines 
(H/I, [Om], and Hu) at a signal-to-noise ratio of ~10 and de¬ 
rive a more accurate measure of the redshift; possibly, this ob¬ 
servation will also infer the presence of kpc-scale outflow in 


the ionized gas component. A more efficient but more time- 
consuming way to detect the outflow on kpc scales would be 
to use VLT/SINFONI in the J and H bands. Moving to sub- 
mm wavelengths, ALMA observations will allow us to extend 
these studies to the molecular component and probe the link be¬ 
tween UFOs and larg e-scale outflows (e.g.. ITombesi et al. 120151: 
iFeruglio et al. Il20l'^ . 
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